Interference Signatures of Abelian and Non-Abelian Aharonov-Bohm effect on 

Neutral Atoms in Optical Lattices 



Ming-Xia Huo, 1 Nie Wei, 1 David A.W. Hutchinson, 1 ' 2 and Leong Chuan Kwek 1 ' 3 

1 Centre for Quantum Technologies, National University of Singapore, 3 Science Drive 2, Singapore 117543 
2 Centre for Quantum Technology, Department of Physics, University of Otago, Dunedin, New Zealand 
3 National Institute of Education and Institute of Advanced Studies, 
Nanyang Technological University, 1 Nanyang Walk, Singapore 637616 

We propose a scheme to generate an effective Abelian U(l) or non-Abelian SU(2) gauge field for 
cold neutral atoms in a ring- or square-shaped optical lattice by using Laguerre- Gauss lasers. The 
synthetic field produced is strongly localized allowing us to study the Aharonov-Bohm effect on the 
neutral atoms. By preparing a coherent state of atoms initially and allowing them to evolve along 
two different paths enclosing the generated magnetic field, we obtain interference signatures of the 
Aharonov-Bohm effect with distinctly different patterns in the detection area for systems exposed 
to a zero, an Abelian U(l) or a non-Abelian SU(2) gauge field. 
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The vector potential, introduced by Maxwell and 
Thomson as a mathematical tool to describe Faraday's 
field, changed ontologically with the discovery of the 
Aharonov-Bohm (AB) effect [1 . The AB phase is the 
archetypical example of a geometric phase as introduced 
by Berry [2]. As such it has spawn many experimental 
and theoretical studies and remains of intense interest. 
Numerous realizations of the AB effect have been studied 
[3]-[5]. Furthermore, in condensed matter, several AB-like 
effects for charged particles have been investigated [6] , in- 
cluding graphene nanostructures [7], quantum rings, dots 
and mesoscopic devices [8j [9] , and carbon nanotubes [10] . 

Ultracold atoms, with their unprecedented control over 
both internal and external degrees of freedom, in optical 
lattices offer a further opportunity to investigate these 
phenomena by utilizing synthetic gauge fields. These 
synthetic fields can be produced in various ways includ- 
ing rotation of the confining trap [TTJ [12] , using adiabatic 
motion of multilevel atoms in a dark-state picture [15] . 
employing spatially dependent two-photon dressing [14], 
and using staggered-aligned states [15] or rotations [16] 
for lattice gauge theories. Further additional terms such 
as spin-orbit coupling have also been added to the emu- 
lated Hamiltonian [T7] . 

Motivated by these studies, we propose a scheme to 
generate Abelian and non-Abelian fields, utilizing laser- 
assisted tunneling [18 . This scheme generates an almost 
perfect solenoid and hence we can investigate theoreti- 
cally the AB effect and its non-Abelian analogue in the 
presence of these fields. In addition, we note that non- 
Abelian gauge fields have been mooted as possible sources 
of topological order of relevance in the context of topo- 
logical insulators p~9j [20] and topological quantum com- 
puting [2T] , 

Model Setup and Basic Principle: We propose a scheme 
to employ Laguerre-Gauss (LG) laser beams to assist 
inter-site tunneling and realize an artificial gauge poten- 




FIG. 1: Schematic diagram for generating effective magnetic 
fields with cold atoms trapped in square (a) and ring (b) lat- 
tices. The atoms can hop from one site to a neighboring 
site via laser-assisted tunneling by shinning LG lasers per- 
pendicular to the lattice surface (a). The LG beams drive the 
spin-conserved transitions between two identical sublevels in 
neighboring sites in (d) to generate Abelian fields, and spin- 
flipping transitions between two different sublevels in neigh- 
boring sites in (e) to generate SU(2) non-Abelian fields. The 
resulting magnetic flux as shown in (c) is non-zero when going 
around a loop including the laser center which is addressed in 
the lattice center and zero when the laser center is excluded 
from the loop. 

tial in both a square lattice Fig. [IJa) and in a ring geom- 
etry Fig. []Jb) [22 . A natural property of the LG laser is 
a spatially dependent orbital angular momentum along a 
closed loop. Our scheme therefore provides an extremely 
simple alternative way to produce a non-vanishing ef- 
fective magnetic, or more exotic, gauge field for neutral 
atoms, without the need for rotation, lattice accelera- 
tions, or the stringent external field conditions necessary 
to utilize degenerate dark states. Apart from the usual 
square lattice as shown in Fig. [l](a) which can be created 
by the interference of counter-propagating laser beams, 
a ring lattice in (b) can also be created by interfering an 
off-resonant LG laser and a plane wave E^e 1 ^ 01-1 * ^. In 
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cylindrical coordinates, an LG laser mode with frequency 
c^o, wave vector fc > an d amplitude Eq propagating along 
the z axis can be denoted by 
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where f pl (r) = (-if x ^MZ^L^e 

r w is the waist of the beam, and Lp' are the Laguerre 
functions [22 . The labels p and I represent the radial and 
azimuthal quantum numbers, respectively. The interfer- 
ence of the LG beam and the plane wave produce the 
periodic interferogram in cpo with Iq wells. After loading 
atoms into the square lattice in (a) or the ID ring in (b), 
we apply Raman LG lasers Ef p i (r) e ll(p e^t-kz) p er p en _ 
dicular to the lattice, as illustrated in (a), to drive atomic 
hopping from site to site. We note that any laser-assisted 
tunneling scheme can be used for this step. For exam- 
ple, one can trap atoms in a hyp erfine- level- independent 
superlattice V (~r^) = Vo J2 a [ cos2 ( nr a) + c °s 2 (27rr a )] 
with ~ 50 — 150kHz. Using 40 K atoms with hyper- 
fine levels F — | as effective spins in the zero-energy 
main sites and F = | as intermediate links with energy 
around 50 — 100kHz [23 , the coherent LG lasers induce 
tunneling of spins F = | along a ring as shown in Fig. 
fl[c). Utilizing one Zeeman level in F = | generates an 
Abelian field [figure 1(d)] whilst using two Zeeman lev- 
els generates SU(2) non-Abelian field [figure 1(e)]. We 
can also use the setups in [15], where two-electron atoms 
generally have a spin-singlet ground state g and a long- 
lived spin triplet excited state e with opposite polariz- 
abilities at some "anti-magic" wavelength, experiencing 
potentials with opposite signs and therefore aligning in a 
staggered manner. After loading, the resonant LG lasers 
then drive tunnelings between staggered states along a 
ring asg—^e—^g—?>e— In this case, considering 
two inner states in g and e also generates the SU(2) field 
as shown in Fig. []Je). 

SU(N) Gauge Field: In lattice Abelian gauge theories, 
the gauge transformation is defined as an U(l) trans- 
formation of the field operators if 0® —> Oi = e lXi O® 

and the phase <^ • acquired in tunneling e^.j (0^0® 
in the Hamiltonian is modified to <fiij — (j)^ - + — 
Xj. Under such gauge fields, the unitary transforma- 
tion describing a particle hopping around a LG driv- 
ing laser is given by the cumulative phase <p = </>i,2 + 
02,3 + 03,4 + • • • 0iv L -i,iV L + 0a/l,i> or the Wilson loop 
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the number of lattice sites. This cumulative phase is re- 
lated to the magnetic flux through the loop via the AB ef- 
fect. When the center of the Raman laser is addressed to 
the middle of a plaquette or the center of the ring, numer- 
ical simulations for square- and ring-shaped lattices show 
that a non-uniform flux with accumulated phase tt is gen- 
erated for I = 1, in the form of an effective solenoid. The 
flux of the solenoid is entirely located within the plaque- 



tte or ring. We further note that even if the center of the 
laser is slightly shifted from that of the lattice, it will not 
affect the results in the thermodynamic limit, although 
some small fluctuations may appear around the solenoid. 
The phase generated in this manner corresponds to a vec- 
tor potential of the form = (Aa.,£c> Aa, v ,0). With the 
lattice sites represented by (m, n) and lattice constant a, 
the vector potential can be written as 
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in units of — according to the standard lattice 
gauge theory [24], where c is the light speed 
and e is an effective charge. Eq. (TTJ de- 
fines the mean Berry connection [TT]. (fm, n = 
/ wl (~r* - f^m+i.n) <P ( x i V) w i - ~^m,n) d 2 ^ and 
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fw^(~r^- r^ m , n+ i) <p (x, y) wi (7^ - r^ m , n ) d 2 ^ 
with cp (x, y) representing the phase of the Raman LG 
beam in the point (x,y). 

When extending the SU (1) transformation to SU (N) 
by considering N inner states, the operator transforma- 
tion becomes 0° 0* = 7*0? with T { a N x N matrix, 
and the Wilson loop W giving the state transformation 
as |^) W |^) becomes W W lj2 + W 2 ,3 + W 3>4 + . . . + 
WOv l -i,n l +Wjnt l ,i with Wij = TjwVjTj. When the loop 
operator is reduced to a global phase W = e %e I with / 
the N x N unit matrix, or equivalently the magnitude of 
Wilson loop is AT, the system reduces to an Abelian sys- 
tem with N independent evolution equations for the N 
"flavor" components. For generating SU(2) fields, as for 
the example shown in Fig. []Je), two degeneracies within 
the Zeeman manifolds should be considered. For exam- 
ple, for 40 K atoms, mp = §, \ in F = § and mp = f , § 
in F = | represent two "colors" of the gauge field, and 
which become split when a moderate magnetic field is 
applied. Here we use a left-circularly polarized LG mode 
with p = 0, I = 1, and a superposition of two right- 
circularly polarized LG modes with p = 0, I = and 
p = 1, I = together to drive the spin-flipping transi- 
tions. We choose the beam width as w/y/2 = N^a/ (2y/2) 
for I = laser and w = for I = 1 laser, such that 

the amplitudes f (r) = \J\e^{^) 2 and /i(r) = (^) 
of lasers with angular momentums I — 0, 1 are approxi- 
mately equal, which gives a unitary hopping matrix and 
generate a non-Abelian SU(2) field. The corresponding 
vector potential ^na = (^-na,^, ^.na,^ 0) is given by 
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The unitary tunneling operators are 
U (m, n — » m + 1, n) = U x (n), U (m + 1, n — >• m, n) = 
JTj (n), C7 (m, n — >• m, n + 1) = U y (m) and 
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FIG. 2: The time evolution of particle distributions for par- 
ticles hopping around a loop formed by LG beams in a ring 
with iVx = 100 sites under a zero gauge field (a), an Abelian 
U(l) field (b), or a non-Abelian SU(2) field ((c) and (d)). 
The particle densities are represented by the percentages of 
particle numbers to the total number of the initially prepared 
Bose-Einstein condensate. At time t — 0, the particles are 
prepared around the site x = 1. During the time period from 
t = to t — T — 600 J -1 , the fringe patterns around the 
opposite site x = TVl/2 become completely different for the 
three field cases (B = 0, U(l) and SU(2)). Around this site, 
there is always a constructive interference for the zero field 
case in (a) and a destructive interference in (b). In contrast, 
for the non-Abelian case, the charge wave density, the sum 
of the two effective spins, has no unique behavior as shown 
in (c), and the spin wave density, as a difference of the two 
effective spins, exhibits two kinds of spin waves as shown in 
(d), where the red-color (blue-color) spin wave has a construc- 
tive (destructive) interference around the site x = and 
a destructive (constructive) interference around the nearest 
neighboring area. 



m) 



where U x 



U (m, n + 1 — >• ra, n) - 
exp (— ieaA-NA,x/fc c ) an d U y = exp (— ieaA^^^/hc). 
The appearance of the phase factor in 
U x (n) U y (ra + 1) Ul (n + 1) (ra) can be under- 
stood in terms of the AB phase accumulated by the 
particles when they travel along a closed contour [TT] . 

Dynamical evolution and AB effect: We load a Bose- 
Einstein condensate (BEC) around one site, away from 
the center of LG-laser, and allow the system to evolve 
under different gauge fields. We then measure the fringe 
patterns with time. For a ID ring-shaped optical lattice 
as illustrated in Fig. [ljb) , numerical simulations for spa- 
tial distributions of particle numbers are plotted in Fig. 
[2j where the particle numbers are shown as a percent- 
age to the total number in the initially prepared BEC. 
We have chosen the lattice site number as = 100 
and the time period T as 600 in unit of J -1 , where the 
typical hopping J ~ 0.05i^R with Er the recoil energy. 



Distinctly different interference pattens can be seen for 
particles experiencing a zero magnetic field as shown in 

(a) , an Abelian U(l) gauge field as shown in (b), and a 
non-Abelian SU(2) gauge field as shown in (c)-(d). The 
reason for the two subfigures in the SU(2) case is that 
two effective spins are involved in the SU(2) field. We 
plot an effective charge wave density (the sum of den- 
sities of the two effective spins) in (c) and an effective 
spin wave density (the difference of densities of the two 
effective spins) in (d), respectively. 

At time t = 0, a BEC is prepared around site x = 
1 = TVl + 1, as shown in Fig. [2] After the initial state 
is prepared, we turn on the LG beams. For the zero- 
magnetic field case in Fig. [2|a), we turn on a laser with 
I = 0, and for the U(l) Abelian field case in (b), we 
turn on a I = 1 laser as illustrated in Fig. [TJd), and for 
the non-Abelian SU(2) field case in (c) and (d), we turn 
on two orthogonally polarized I = and I = 1 lasers as 
illustrated in Fig. [I] (e). After a certain time, we see 
that the particle number distributions around the site 
x = directly opposite to the initially occupied site, 

are different in the three cases [see Fig. [2](a)-(c)]. Note 
that a destructive interference is always seen in Fig. [2] 

(b) at site The spin wave density in (d) can be 
regarded as a difference of two effective spin densities: 
one colored red and the other blue in Fig. [5] (d). At 
site the constructive interference is always due to 
one type (red). This result can also be explained as a 
superposition of two effective spin waves under a zero 
gauge field due to the I = laser and a finite gauge field 
due to the I = 1 laser. The red-color spin wave therefore 
represents the spin wave sensitive to the zero magnetic 
field and the blue-color spin wave represents the other 
spin wave sensitive to the non-zero magnetic field. 

Our key result for the ring geometry is a clear demon- 
stration of the AB effect. In the presence of an Abelian 
synthetic gauge field, the AB phase always induces de- 
structive interference at site Absence of any popu- 
lation at site A/l/2 at any time is therefore a clear signa- 
ture of the Abelian AB effect. Furthermore, we demon- 
strate the effects of the non-Abelian gauge field in terms 
of the spin density waves. 

For the 2D square lattice case as illustrated in Fig. [TJa) 
andJTJc), we can also prepare the initial state by load- 
ing a BEC around one plaquette. For the U(l) Abelian 
field case as shown in Fig. |3jb), we turn on a I = 1 
laser as illustrated in Fig. [ljd) , and for the non-Abelian 
SU(2) field case as shown in Fig. [3] (c)-(d), we turn on 
two orthogonally polarized I = and I — 1 lasers as il- 
lustrated in Fig. [lje). Numerical simulations with time 
show that the angular momentum imparted by the LG 
beams results in a circular distributions of atoms. The 
particle number distribution at the opposite side of the 
circle from the initially prepared site again shows clear 
destructive interference in the case of the Abelian gauge 
field - the signature of the AB effect. As a concrete ex- 
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FIG. 3: The time evolution of particle distributions for par- 
ticles hopping around a loop formed by LG beams in a 2D 
square lattice with 40 x 40 sites. At time t = as shown 
in (a), the particles are distributed around one site. The 
screenshots at time t = 2000 J -1 are given for the system 
under an Abelian U(l) (b) or a non-Abelian SU(2) ((c) and 
(d)) gauge field. The particle densities are represented by the 
percentages of particle numbers to the total number of the 
initially prepared Bose-Einstein condensate. While a destruc- 
tive interference always occurs around the opposite site in the 
Abelian field case (b) , a non- vanishing charge wave density oc- 
curs everywhere in (c) and two kinds of spin waves appear in 
(d) as red- and blue-colored spin waves, where the red-color 
(blue-color) spin wave has a constructive (destructive) inter- 
ference at the opposite side and a destructive (constructive) 
interference at its two neighboring area. 



ample, screenshots at time t = 4000 J 1 are given in Fig. 
[3] (b)-(d). The effect is similar to the ring geometry case. 

In summary, we have presented a method for generat- 
ing Abelian and non-Abelian synthetic gauge fields with 
cold atoms. As an example of the use of this method, 
we demonstrate an interferometric scheme to observe di- 
rectly the AB effect. We also illustrate interference pat- 
terns for charge and spin density waves corresponding 
to the non-Abelian AB analogue. These effects are ro- 
bust against decoherence as they are topological in ori- 
gin. Consequently, they could in principle be harnessed 
for quantum information processing [25 -27 . We have 
demonstrated a method for generating U(l) and SU(2) 
gauge fields, this technique in principle can be general- 
ized to SU(JV). 

We acknowledge support from the National Research 
Foundation & Ministry of Education, Singapore. 
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